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Combating the effects of disorder on light transport in micro- and nano-integrated 
photonic devices is of major importance from both fundamental and applied view¬ 
points. In ordinary waveguides, imperfections and disorder cause unwanted back- 
reflections, which hinder large-scale optical integration. Topological photonic struc¬ 
tures, a new class of optical systems inspired by quantum Hall effect and topological 
insulators, can realize robust transport via topologically-protected unidirectional edge 
modes. Such waveguides are realized by the introduction of synthetic gauge fields for 
photons in a two-dimensional structure, which break time reversal symmetry and en¬ 
able one-way guiding at the edge of the medium. Here we suggest a different route 
toward robust transport of light in lower-dimensional (ID) photonic lattices, in which 
time reversal symmetry is broken because of the non-Hermitian nature of transport. 

While a forward propagating mode in the lattice is amplified, the corresponding back¬ 
ward propagating mode is damped, thus resulting in an asymmetric transport that is 
rather insensitive to disorder or imperfections in the structure. Non-Hermitian trans¬ 
port in two lattice models is considered: a tight-binding lattice with an imaginary 
gauge field (Hatano-Nelson model), and a non-Hermitian driven binary lattice. In the 
former case transport in spite of disorder is ensured by a mobility edge that arises 
because of a non-Hermitian delocalization transition. The possibility to observe non- 
Hermitian delocalization induced by a synthetic ’imaginary’ gauge field is suggested 
using an engineered coupled-resonator optical waveguide (CROW) structure. 
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The realization of photonic devices at the micro- and nano-scale capable of guiding light in a controllable way 
despite the presence of disorder or imperfections is of major importance for integrated optic applications. Inspired 
by the quantum Hall effect and the concept of topological insulators in condensed-matter physics, several recent 
works have proposed and experimentally demonstrated one-way edge propagation for light waves in two-dimensional 
photonic structures (for a recent review see [I]). Topologically protected edge states are one-way guided modes 
propagating (clockwise or counter-clockwise) along the edge of the sample that cannot be scattered into other states, 
and are therefore immune to back reflection and localization. The existence of such states requires to break time 
reversal symmetry of the optical system using magnetic materials, such as in magnetic photonic crystal systems [Ml, 
or by creating synthetic gauge fields for photons. The latter have been demonstrated in a wide variety of optical 
systems, such as coupled silicon resonators 15 3 - twisted waveguide lattices 0, bianisotropic nretamaterials [9], and 
lattices with parameters modulated in time m- The possibility to observe and adiabatically transport edge states in 
quasi-crystals has been reported as well In¬ 
in this work we suggest a different route toward robust propagation in one-dimensional (ID) photonic lattices, 
in which time reversal symmetry breaking is realized by exploiting non-Hermitian dynamics. Wave propagation in 
non-Hermitian (complex) crystals can show a wide variety of interesting effects as compared to ordinary (Hermitian) 
crystals, such as double refraction and nonreciprocal diffraction [12] , unidirectional Bragg scattering and invisibility 
|l3lfl5| . anomalous dynamic localization |16| . hyperballistic transport E2. and a transition from ballistic to diffusive 
regimes [18]. Such effects have been mainly investigated in case of VT symmetric systems pjj. However, in such 
systems wave transmission turns out to be independent of the incidence side, so VT structures are unsuited to 
realize robust one-way transport. Indeed, asymmetric transmission and optical isolation here requires to exploit e.g. 
nonlinearity in the system [[20] l21j . Other possibilities to introduce asymmetric transport is the use of ’imaginary’ 
gauge fields or modulation in time of some parameters of the optical medium. In a pioneering paper, Hatano and 
Nelson [22] showed that an ’imaginary’ magnetic field in a disordered ID lattice can induce a delocalization transition, 
i.e. it can prevent Anderson localization (see also E3H13I13)- In the disordered lattice, delocalized states correspond 
to complex energies, whereas real energies are associated with Anderson-like localized states. This is a very interesting 
result because it shows that mobility is not prevented by disorder in the non-Hermitian ID lattice with an ’imaginary’ 
gauge field, contrary to a ’real’ magnetic held that has no effect on Anderson localization in a ID chain. Such a result, 
however, has been overlooked in the condensed matter physics context, because of the challenging task of synthesizing 
an ’imaginary’ vector potential. In optics, however, this can be realized in coupled microring structures, as we will 
discuss in our work. We also suggest a more general and physically simple understanding of robust transport in 
non-Hermitian lattices with asymmetric transmission, the imaginary gauge potential being a special example. 


RESULTS 

One-way transport in non-Hermitian photonic lattices: Introductory physical idea. One-way propagation 
in 2D Hermitian lattices with a synthetic gauge held is based on the existence of topologically-protected edge states. 
Here we suggest a different mechanism of one-way robust propagation in a non-Hermitian ID lattice, in which a 
preferred propagation direction arises because of unbalanced amplihcation/attenuation of waves propagating in op¬ 
posite directions. The main physical distinction between topologically protection of edge waveguides in a Hermitian 
2D lattice and asymmetric transport in a non-Hermitian lattice can be readily understood in reciprocal space, as 
schematically illustrated in Fig.l. In an Hermitian 2D lattice with a synthetic gauge held, the dispersion diagram 
E w = E w {q ) of the edge guided mode shows the typical behavior depicted in Fig. 1(a). It connects the bulk frequency 
bands above and below the frequency gap. The unidirectionality of the edge guided mode can be seen from the 
group velocity v g = ( dE w /dq ) of the waveguide dispersion curve, which has only positive (or only negative) value. In 
addition, there are no counter-propagating modes at the same frequencies as the one-way edge mode. This enables 
light to how around imperfections with perfect transmission. Let us now consider a non-Hermitian ID lattice with a 
transmission band described by the energy dispersion curve E = E(q), such that the following condition is satished 

v g (q)lm.(E{q)) > 0 (1) 

where q is the Bloch wave number, which varies in the interval — tt < q < 7r, and v g (q) = Re(dE/dq) is the group 
velocity [Fig. 1(b)]. In this case, a forward propagating mode, corresponding to v g > 0, is amplified, whereas the 
corresponding backward propagating mode, with a group velocity v g < 0, is damped. This unbalanced amplihca¬ 
tion/attenuation of waves introduces a preferred direction of propagation, and back scattered waves arising from 
disorder or imperfections in the lattice experience attenuation. Thus localization effects are expected to be prevented. 

Imaginary gauge field and non-Hermitian delocalization. The simplest example of a non-Hermitian ID tight- 
binding lattice with a dispersion relation satisfying the condition (1) is provided by the Hatano-Nelson model [25] . 
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Originally introduced to study the motion of magnetic flux lines in disordered type-II superconductors, it describes 
the hopping dynamics of a quantum particle on a tight-binding ring lattice threaded by an imaginary magnetic flux. 
In the absence of disorder, the corresponding Hamiltonian reads [53] EJ] 

H = ^ k {exp(— h)\n)(n + 1| + exp(/i)|n + l)(n|} (2) 

n 

where k is the hopping rate between adjacent sites and the parameter h describes the effect of an imaginary vectorial 
potential. A simple physical implementation of an imaginary gauge field for photons in a CROW structure will be 
discussed in details in a following section. Here we note that, for an infinitely-extended lattice, the dispersion relation 
of the tight-binding lattice band is given by 

E(q) = 2 k cos (q + ih) = 2 k\ cos q — 2zk- 2 sin q (3) 

where we have set K\ = Kcosh h and K 2 = k sinh h. The group velocity is given by v g = R e(dE/dq) = — 2Kisin<7. 
Note that a forward-propagating wave (—7r < q < 0, v g > 0) is amplified because lm(E(q)) > 0, whereas a backward 
propagating wave (0 < q < 7r, v g < 0) is attenuated because Im(E(q)) < 0. This makes wave transport in the lattice 
highly asymmetric because backward propagating waves vanish after some propagation distance due to damping. As 
a result, unidirectional transport is expected to be robust against disorder and lattice imperfections. For disorder in 
the site energies, the Hamiltonian takes the form 

H = ^ k {exp(—h)\n)(n + 1| + exp(/i)|n + l)(n|} + ^2Vn\n)(n\ (4) 

n n 

where V n is a random potential. V n is assumed independently for each site, from a distribution with zero mean and 
variance u 2 . After setting \il>(t)) = 'Yh n c n {i)\n), the evolution equations for the amplitude probabilities c n (t) read 
explicitly 


,dc r 

i—r 

dt 


Kexp(-h)c n+1 + Kexp(/i)c„_i + V n c n . 


(5) 


For h = 0, it is known that all the eigenstates of H become localized owing to Anderson localization, regardless of 
the strength u 2 of disorder. However, for ft / 0 some delocalized states survive, i.e. a delocalization transition is 
observed [551155] . Rather generally, the delocalized (localized) states correspond to complex (real) energies, however 
the precise nature of the non-Hermitian eigenstates is somewhat controversial [241 ; 55] . As discussed in previous 
works, by introducing an imaginary vector potential a mobility region opens up near the center of the tight-biding 
energy band, where Anderson localized states show the largest localization length, i.e. they are weakly localized by 
disorder. Such a phenomenon can be explained as follows [55] . The eigenstates of the Anderson lattice with the 
imaginary vector potential can be mapped into the ones of the Hermitian Anderson lattice (i.e. with h = 0) after 
the transformation c n = a n exp(hn). Hence, if a n is a localized state in the Hermitian Anderson lattice with energy 
E = Eq and localization length £ = £(£))), it remains a localized state in the non-Hermitian Anderson lattice provided 
that 2\h\ < l/£(-E 0 ). Such a condition in violated in the range of energies —E c < E 0 < E c around E 0 = 0, i.e. a 
mobility region is opened with an edge at E c . The mobility edge E c can be computed following the analysis of Ref. 
[23] and reads explicitly E c ~ (4k 2 — u 2 /2\h\) 1 / 2 . The mobility region opened by the imaginary gauge field ensure 
particle propagation along the lattice in spite of disorder. This is shown, as an example, in Fig.2. In Figs.2(a) and 
(b) a semi-infinite lattice is considered, and the edge site at n = 0 is occupied at initial time t = 0. The figure shows 
the numerically-computed evolution of the wave packet center of mass 




E n n\c n (t )\ 2 

Zn\Cn(t )\ 2 
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for the Hermitian [Fig.2(a)] and non-Hermitian [Fig.2(b)] lattices. As in the former case ( n{t)) clearly ceases to 
increase as t increases (this is a clear signature of Anderson localization), in the latter case a secular growth of ( n(t )) 
is observed, indicating the existence of delocalized (mobility) states. The robustness of the transport against disorder 
in the non-Hermitian lattice is also observed when the lattice is initially excited in a single site in the bulk, i.e. far 
from the lattice boundaries, or by a localized wave packet. As an example, Fig.2(c) showa the numerically-computed 
dynamics in the disordered lattice with initial excitation of site n = 100. Since the non-Hermitian dynamics does 
not conserve the norm, for the sake of clearness in the figure the evolution of the site occupation probabilities for 
the normalized amplitude probabilities a n (t) = c n (t)/\c n (t)\ 2 is depicted. The figure clearly shows that, in 
spite of disorder, transport is observed in the forward direction. The same scenario is found for initial wave packet 
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excitations, as discussed in the Methods. In this case, even in the absence of disorder, owing to the dependence of 
the complex energy dispersion curve E(q) on the Bloch wave number q [Eq.(3)], during the propagation an initial 
wave packet of the form c n (0) = F(n) exp(iq 0 n), with carrier Bloch wave number q 0 and slowly-varying amplitude 
F(n), suffers for reshaping (distortion) effects because of both group velocity (phase) and amplification (amplitude) 
dispersion, i.e. because R e(d 2 E/dq 2 ) / 0 and Im (dE/dq) ^ 0. For the Hatano-Nelson model, both group velocity and 
amplitude dispersion effects are minimized at qo = 7t/ 2, where R e(d 2 E/dq 2 ) = lva.{dE/dq) = 0. Examples of wave 
packet propagation are discussed in the Methods. Finally, it should be pointed out that non-Hermitian transport is 
robust also against structural imperfections or defects in the lattice. Let us consider, as an example, a lattice with 
two potential defects at sites no and ni, i.e. let us assume in Eq.(4) V n = Vo {S n ,n 0 +^n,ni)i where Vo is the strength of 
the potential defect. In the Hermitian lattice (h = 0), a propagating wave packet undergoes multiple reflections back 
and forth between the two defects, like in a Fabry-Perot cavity. This yields multiple transmitted wave packets, i.e. 
echoes of the original wave packet, as illustrated in Fig.3(a). Application of the imaginary gauge field to the lattice 
(h ^ 0) suppresses multiple reflections (echo effects), as shown in Fig.3(b). 

Non-Hermitian driven lattice model. We conjecture that the absence of Anderson localization and asymmetric 
transport induced by an imaginary gauge field is a general feature of a non-Hermitian ID lattice with a dispersion 
relation of the energy band that satisfies Eq.(l), i.e. that provides amplification (attenuation) for forward (backward) 
propagating waves. In a ID lattice with local complex site potentials, transport is always symmetric (reciprocal) 
126] . However, non-reciprocal (asymmetric) transport can be obtained by introduction of some periodic modulation 
of system parameters. In this case, the quasi-energy bands of the driven lattice may satisfy Eq.(l), i.e. asymmetric 
transmission can arise. To support our conjecture, let us consider as an example the periodically-driven non-Hermitian 
binary lattice described by the Hamiltonian 

h = K i\ n )( n + !| + \ n + iX"!} + { Fn + (-irvK*)} \ n )( n \ + v n\ n )( n \ ( 7 ) 

n n n 

where n is the hopping rate between adjacent sites, F is an index gradient along the lattice, V n is a random potential 
that accounts for disorder, and 


g(t) = Gr cos (cot) + iGi sin(wf) (8) 

is a sinusoidal function of frequency u> that describes a periodic modulation of the complex energy sites with opposite 
sign at alternating sites. Such a non-Hermitian binary lattice can describe temporal hopping of light in a CROW 
structure with an impressed static index gradient and with temporal modulation of gain/loss at altering sites in the 
chain JSJ. Also, provided that the temporal dynamics is replaced by spatial propagation, Eq.(7) can describe spatial 
light transport in an array of circularly-curved evanescently-coupled optical waveguides with alternating gain/loss and 
index modulation in the guides along the propagation direction t nzi!2g]. The disorder V n accounts for imperfections 
in the propagation constants of the various waveguides in the array. Note that the lattice becomes Hermitian in the 
Gj = 0 limit. In the absence of disorder, H is periodic in time with period T = 2tt/u, and Floquet theory applies. 
Provided that the resonance condition F = ( M/N)u> is satisfied for some irreducible integers M and N, the lattice 
sustains two quasi-energy minibands with dispersion relations that satisfy Eq.(l) (see Methods). As an example, 
Fig.4(a) shows the numerically-computed quasi energy spectrum of the two mini bands for ui/k = 1, F/w = 2, 
Gr/k = 4.7 and Gi/k = 4.26. The quasi-energy spectrum satisfies the condition (1) and hence forward (backward) 
waves are amplified (damped) as they propagate in the modulated lattice. The robustness of forward-propagation in 
the presence of disorder is demonstrated in Figs.4(b) and (c). A semi-infinite lattice in considered, and excitation of 
the edge site is accomplished at initial time t = 0. Figure 4(b) depicts the evolution of the wave packet center of mass 
(n(t)) for the driven lattice in the Hermitian limit Gi = 0, clearly showing localization in the presence of disorder. 
Conversely, in the non-Hermitian case propagation is not prevented by the disorder [Fig.4(c)]. A similar scenario 
is found by considering initial excitation of a site in the bulk. Robust transport occurs for an initial distribution 
corresponding to a localized wave packet as well (see Methods). Interestingly, as compared to the Hatano-Nelson 
model, in the driven lattice model (5) wave packet distortion effects arising from group velocity and amplification 
dispersion are negligible because in this case one has Re(d 2 E / dq 2 ) ~ 0 and lm(dE/dq) ~ 0 almost everywhere inside 
the Brillouin zone [Fig.4(a)]. 

Realization of an imaginary gauge field in coupled resonator optical waveguides and non-Hermitian 
delocalization. Synthetic gauge fields and robust transport via topologically-protected edge modes in 2D CROW 
structures have been proposed and experimentally demonstrated in a series of recent works using coupled microring 
resonators m- Here we show that, by inclusion of engineered gain and loss terms in auxiliary ring resonators, a 
synthetic ’imaginary’ gauge field can be readily implemented in a linear (ID) chain of microring resonators, which 
can thus provide an accessible optical system for the observation of non-Hermitian delocalization [22] . A schematic of 
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the ID chain of coupled ring resonators is shown in Fig.5(a). The main ring resonators are indirectly coupled using 
another set of auxiliary rings which are designed to be antiresonant to the main ring resonators, i.e., the length of 
the connecting rings is slightly larger (or smaller) than the main rings so as to acquire an extra 7 r phase shift. The 
auxiliary ring provides amplification in the upper half perimeter, with single-pass amplification h , and balanced loss in 
the lower half perimeter, with single-pass attenuation —h. Indicating by a n (t) the amplitude of the counterclockwise 
propagating field in the n-th ring in the main resonators, with a carrier frequency coincident with one longitudinal 
ring resonance, coupled-mode equations for the slowly-varying amplitudes a n (t) can be derived in the mean-field limit 
after elimination of the field amplitudes in the auxiliary resonators (see Methods). In the absence of disorder the 
equations read 

T ~^F = _7a ” _ iK { ex P(“ /l ) a «+i + exp(/i)a n _i} (9) 

where r is the round-trip time in the main rings, 7 their single-pass loss rate, and n is an effective coupling rate, 
defined by Eq.(26) given in the Methods. Note that, after the substitution a n (t) = c n (t) exp(—yt/r), Eqs.(9) reduces 
to the Hatano-Nelson model (5) without disorder. It should be noted that, if the circulation of the light fields in the 
rings is reversed, Eq.(9) still holds after the substitution h —► —h, i.e. the direction of robust transport is reversed. 
However, provided that light circulation in either one of the two directions is not excited or prevented, the CROW 
structure realizes the non-reciprocal Hatano-Nelson model (5). In the CROW system, the disorder V n arises because of 
deviations of microring resonance frequencies from the reference value due to fabrication imperfections. Deviations of 
the antiresonance condition in the auxiliary rings also contribute to the disorder V n (see Methods). In the Hermitian 
limit h = 0, i.e. in the absence of the gain/loss regions in the auxiliary rings, disorder is responsible for Anderson 
localization, which has been experimentally observed in Ref. [31] . The introduction of the imaginary gauge field should 
prevent Anderson localization. 

Here we suggest a simple experiment for the observation of the non-Hermitian delocalization transition that occurs 
in the presence of the imaginary magnetic field. To this aim, let us consider a finite CROW chain made of N 1 
main microrings, in which the ring in one edge of the chain (for instance the ring with index n = 1 ) is pumped and 
provides a gain parameter go per round-trip; see Fig. 6 (a). Taking into account disorder of the microring resonance 
frequencies and the gain in the edge microring, coupled-mode equations are modified as follows 

= + ~ iV n )a n - in {exp(—/i)a n +i + exp(/i)a„_i} (10) 

(n = 1,2,3, ...,1V). The CROW structure sustains a set of N modes (also called supermodes), which are obtained 
as the eigenvectors of the linear system (10). For a given realization of disorder V n and for go = 0, all modes are 
damped in time, because their associated eigenvalues have a negative real part. As the gain parameter g 0 is increased 
up to some threshold value, a neutrally stable (i.e. undamped) mode will arise among the N supermodes, which 
corresponds to the lasing mode of the CROW structure. The localization properties of the lasing mode turn out to 
be strongly sensitive on the value of h, i.e. of the imaginary gauge field. In fact, the N supermodes sustained by the 
active microring chain and their corresponding thresholds go are obtained from Eq.(10) after making the Ansatz 

a n (t) =p n exp(zAt/r), ( 11 ) 

where A is a real parameter that defines the frequency detuning of the supermode from the microring reference 
frequency. Substitution of Eq.(ll) into Eq.(10) yields 

(iX + 7 - g 0 Sn,i)Pn = -iVnPn - in {exp(-/i)p n+ i + exp(ft)p n _i} . ( 12 ) 

If the gain (go) and loss ( 7 ) parameters are smaller or of the same order than /c, the mode distribution \p n \' 2 is not much 
distorted from the one of the chain with 7 = g 0 =0. Hence the localization property of the lasing mode is sensitive to 
the absence or the presence of the imaginary gauge field h. For h = 0 and for a sufficiently strong disorder, because 
of Anderson localization all the modes are localized in a few sites of the chain with a localization length £ < N, and 
the lasing mode corresponds to the localized state p n (h = 0) with the largest occupation at the pumped ring n = 1 
[Figs. 6 (b) and (c)]. For h / 0, the distribution of the lasing mode is simply given by p n (h) = p n (h = 0) exp (/in), as 
it readily follows from Eq.(12) after the transformation p n —> p n exp(/m). Hence, as h is increased such that 2 h gets 
comparable to the inverse of the localization length l/£, the mode p n (h) is delocalized all along the chain or even 
localized closest to the opposite edge n = N if h is further increased. Hence the oscillating mode will be delocalized 
[Fig. 6 (d)]. In an experiment, the delocalization effect induced by the imaginary gauge field h can be thus simply 
visualized as a delocalization transition of the lasing intensity distribution in the microring chain as the parameter h 
is varied. Such a simple experiment might provide the first accessible testbed for the observation of the non-Hermitian 
delocalization transition predicted in the pioneering work by Hatano and Nelson. 
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DISCUSSION 

Non-Hermitian photonic lattices with asymmetric amplification/attenuation of counter-propagating modes can 
provide a route toward robust light transport, preventing Anderson localization in the presence of disorder. The 
physical origin of one-way transport in non-Hermitian lattices considered in our work is rooted in the non-Hermitian 
delocalization transition originally discussed in quantum mechanics by Hatano and Nelson for a lattice with an 
imaginary gauge potential [22] . It is thus very distinct than the physics of topologically-protected edge modes arising 
in 2D photonic lattices with synthetic gauge fields [I], or the physics of asymmetric transport in other VT -symmetric 
systems where non-reciprocity is obtained by exploiting some nonlinearities in the system mm . We have suggested 
a very simple physical condition for non-Hermitian robust one-way transport and non-Hermitian delocalization [see 
Eq.(l)], and shown that the Hatano Nelson model provides a special example where such a condition is met. We 
conjectured the generality of such a condition, that can be satisfied in driven non-Hermitian lattices as well (see 
Fig.4). Finally, we have proposed a simple method for the implementation of an imaginary gauge field for photons in 
coupled microring resonators, which could thus provide the first experimentally accessible testbed for the observation 
of the Hatano-Nelson non-Hermitian delocalization. Our results indicate that non-Hermitian photonic lattices can 
provide a new route toward unidirectional light transport which is robust to imperfections and disorder, and could 
be of relevance to other related phenomena in photonics, such as non reciprocal photonic transmission, isolation, 
and unidirectional coherent perfect absorption fTH?l-i.‘S5j . as well as in other non-Hermitian physical systems such as 
Josephson circulators [34l 135] . 


METHODS 

Quasi-energy spectrum of the driven non-Hermitian lattice. The quasi-energy spectrum of the Hamiltonian 
(7) in the absence of disorder (V n = 0) can be calculated as follows. After setting \i/>(t)) = ^2 n a n (t)exp(—iFnt)\n), 
the evolution equations for the amplitude probabilities a n read 

= k {exp(-iFt)a n+ i + exp(iFt)a n -i} + (-1 ) n g(t)a n . (13) 

We look for a solution to Eq.(13) of the form 

a n (t) = ^ ^ exp(-ig?r/2) (14) 

where the upper (lower) row applies to an even (odd) value of n, and where q is the Bloch wave number that varies 
in the range (—7r,7r). Substitution of Eq.(14) into Eq.(13) yields 


dA , , , ,„ , 

i— = 2 kB cos (g/2 + Ft) + g(t)A 

(15) 

dB . , , 

i— = 2kAcos (g/2 + Ft) — g{t)B. 

(16) 


If the resonance condition F = ( M/N)u> is satisfied, Eqs.(15) and (16) are periodic in time with period T = 2-kN/uj 
and Floquet theory applies. Indicating by $(g) = exp(— iTTZ(q)) the propagator of Eqs.(15) and (16) from t = 0 to 
t = T, i.e. (A(T), B(T) t = *I>(g)(A(0), B(0)) T , the quasi energies E±(q) of the two mini bands for the Hamiltonian 
H are defined as the eigenvalues of the 2x2 matrix 7 Z(q). Quasi energies are defined apart from integer multiples 
than uj/N; for the sake of definiteness the real part of the quasi energies is taken within the interval (—u>/2N,uj/2N). 
Note that, since the modulation function g(t) satisfies the condition g(—t) = g*(t ), then it can be readily shown that 

E±(-q) = E* ± (q). 

Wave packet distortion effects. Let us consider an initially localized wave packet with carrier wave number 
go that propagates in the ordered lattice. The initial excitation can be written rather generally as c„(0) = c(n,0) = 
F{n) exp(igon) = J dqS(q) exp(ign) with a slowly-varying amplitude F(n ), where S(q) is the Bloch spectrum which 
is assumed to be a narrow function at around q = qo- For a static lattice with a single tight-binding band, like for the 
Hatano-Nelson model (2), the amplitude probabilities c n (t) evolve according to c(n, t) = f* dqS(q) exp[iqn—i.E(q)t] = 
f* dqS(q) exp[ImE(g)f] exp[ign — ?ReE(g)f], where E(q) is the complex energy dispersion curve of the lattice band. A 
similar expression holds for a time-periodic Hamiltonian [e.g. the model (7)], provided that time t is taken at integer 
multiplies than the modulation period 27r/w and E(q) is replaced by the quasi-energy band. Owing to the dependence 
of the complex energy E(q) on the Bloch wave number g, the wave packet is generally distorted during propagation. 
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However, for a wave packet spectrally narrow at around the carrier wave number q = qo such that Im.(dE/dq) qo = 0 
and Re(d 2 E/dq 2 ) qo = 0, at leading order one has \c{n,t)\ 2 ~ exp(2gt)\c(n — v g t, 0)| 2 , where g = Im.E(c/o)> and 
v g = R e(dE/dq) qo . This means that the wave packet propagates nearly undistorted with a group velocity v g , apart 
from a uniform amplification (g > 0) or attenuation (g < 0). In the Hatano-Nelson Hamiltonian the condition of 
nearly-undistorted propagation is attained at qo = 7r/2, whereas for the driven lattice model (7) nearly undistorted 
propagation is predicted for almost any wave number qo, owing to the flatness of lm.E(q) and the linear dependence of 
Re-E(g) on q inside the Brillouin zone [Fig.3(a)]. As an example, Fig.7(a) shows the numerically-computed evolution 
of a two-humped wave packet F(n) = n exp[— (n/10) 2 ] for the Hatano-Nelson Hamiltonian for carrier wave number 
qo = 0 and qo = 7t/2, in either an ordered and a disordered lattice. The figure clearly shows that, according to the 
previous analysis, even in the absence of disorder strong wave packet distortion is observed for qo = 0, whereas wave 
packet distortion is much weaker for q 0 = tt/2. For the modulated lattice Hamiltonian [Eq.(7)], wave packet distortion 
is almost absent, as shown in Fig.7(b). The reason thereof is the special behavior of dispersion curves of the quasi 
energy minibands of the driven lattice, as discussed in the main text. 

Imaginary gauge field in a chain of microrings: coupled-mode analysis. For the sake of clearness, let us 
consider the case of two high-Q rings in the main chain coupled by an auxialiay ring [Fig.4(b)]. The method can be 
readily extended to an arbitrary number of rings in the main cavity indirectly coupled by auxiliary rings. The fields 
in the two rings of the main chain are assumed to propagate counterclockwise, whereas the field in the auxiliary ring 
propagates clockwise. The carrier frequency of the fields is assumed to be one longitudinal frequency of the rings in 
the main chain. Indicating by a(t), b(t ), e(t), f(t), g(t) and l(t) the field amplitudes in the main and 

auxiliary rings at the locations depicted in Fig.4(b), resonator coupling is described by the following equations [SU] 


/ 


u p\(a\ 
-p* u J \ e ) 



u p 
—p* u 



(17) 


where p is the coupling constant and |u | 2 + |p | 2 = 1. We assume small coupling, i.e. |p| <C 1 and u ~ 1. A real value 
of u is also assumed for the sake of clearness. The field boundary conditions for the two main rings read 

a(t + r) = a!(t) exp(— 7 ) ~ (1 — 7 ) 0 !(t) , b(t + r) = b'(t) exp(— 7 ) ~ (1 — fab'(t) (18) 

where 7 <C 1 is the single-pass loss in each ring and r is the transit time in the ring. Taking into account propagative 
effects in the two half sections of the auxiliary ring one can also write 

g(t + ri) = f{t) exp(/i -(- ifa) , e(t + n) = l(t) exp(—h + ifa) (19) 

where 2 ti ~ r is the propagation time in the auxiliary ring, <f>\p, are the accumulated phases due to field propagation 

in the upper and lower sections of the auxiliary ring, and h is the gain/loss parameter in each of the two sections. 

Since the auxiliary ring is assumed to be in anti-resonance, one has fa + fa = ^ (apart from integer multiplies than 
27 t). The anti-resonance condition is obtained by slightly shortening (or lengthening) the length of the auxiliary 
ring as compared to the rings in the main chain. In the single longitudinal mode regime and assuming that the 
field amplitudes vary slowly over one round trip in the rings (mean-field limit), i.e. a(t + r) ~ a(t) + ( da/dt)r , 
b(t + t) ~ b(t) + ( db/dt)r , etc. from Eqs.(17) and (18) one can write 

da . . db . . . 

r— ~ —(1 — u)a — 7 a + pe , r— ~ —(1 — u)b — 70 + pg. ( 20 ) 

dt dt 

On the other hand, from Eqs.(17-19) one has 

g(t + ri) + u 2 g(t — ti) = up*b{t — n) — p*a(t ) exp {h + i(f>\) (21) 

e(t + T\) + u 2 e(t — -77) = up*a(t — 17) — p*b{t) exp {—h + i<j> 2) ( 22 ) 

where we used the anti-resonance condition fa + 0 2 = 7 r. For |p| <C 1, u ~ 1 at leading order the delay effects in the 
amplitudes g and e can be neglected, thus obtaining 

g(t) ~ [ub(t) - a[t) exp(h + ifa)\ , e(t) ~ [-b(t) exp(-/i + ifa) + ua(t)]. (23) 

Assuming that the two gain/loss sections in the auxiliary rings have the same length, i.e. fa = fa = tt/ 2, one obtains 

* 

g(t) ~ [ub(t) — ia(t ) exp(/i)] , e(t) ~ [—ib(t) exp (—h) + ua(t)] . 


(24) 


Substitution of Eqs.(24) into Eqs.(20) yields 


da 


T 


dt 


(l — u — u\p\ 2 /2) a — ja — in exp {—h)b , 



where we have set 


(l — u — it|p| 2 /2) b — 76 — i«exp(/i)a (25) 



(26) 


Taking into account that u = yjl — \p \ 2 ~ 1 — \p\ 2 /2, at leading order from Eqs.(25) one finally obtains 

da db 

r—= — 7 a — iKexp(—/i )6 , r— = — 70 — in exp(/i)a. (27) 

Uib UjL 

Such equations clearly show that the role of the auxiliary ring is to indirectly couple the two main rings, with 
unbalances hopping rates «exp(— h) and «exp(/i), where h is the gain/loss parameter in the two half-sections of the 
auxiliary ring and n is defined by Eq.(26). 

In the above analysis we assumed exact anti-resonance of the auxiliary resonator, i.e. </>i = (f >2 — 7 t/ 2 . A slight 
deviation of the auxiliary ring from the anti-resonance condition can be taken into account by assuming in the previous 
analysis (f>i = 7r/2 + 07 and (f >2 = tt/2 + 1 V 2 , where u >are additional phase deviations from the ideal condition that 
could arise in fabrication imperfections. Anti-resonance is broken whenever wi + W 2 / 0. For small values of phase 
deviations oj\ and W 2 , after some straightforward calculations it can be shown that at leading order Eq.(27) is modified 
as follows 


da 

T 5 = l 11 ~ 7 


■ U2 ^ a — in exp (—h + i0)b , 


—7 + in 1 + — ) b — *«exp(/i — id)a, (28) 


where we have set 8 = ( 0 J 2 — uq)/2. Equation (28) shows that deviation from the anti-resonance condition has two 
effects: it introduces an effective shift of the resonance frequency of the main rings by the amount n(u> \ + u>2)/(2t), 
and it provides a phase term 8 for the hoping rate k. The latter does not play any role since it can be removed by a 
proper gauge transformation. Therefore, slight deviations from the anti-resonance condition of the auxiliary ring can 
be accounted for by a disorder V n of the main ring resonance frequencies, like in Eq.(10). 
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FIGURE CAPTIONS 


Fig.l Principle of robust light transport based on (a) topologically-protected edge states, and (b) asymmetric non-Hermitian 
transport. 

Fig. 2 Numerically-computed evolution of the wave packet center of mass ( n(t )) in a semi-infinite Hermitian lattice de¬ 
scribed by the Hamiltonian (4) corresponding to single-site excitation, (a) Edge site excitation [c n (0) = <5 n ,o] in the Hermitian 
lattice (h = 0), (b) edge-site excitation [c n (0) = <5n,o] of the the Non-Hermitian lattice (h = 0.2), (c) site excitation in the bulk 
\Cn (0) = S„ , ioo] of the the Non-Hermitian lattice (h = 0.2). The dashed curves show the center of mass path in the absence 
of disorder, whereas the solid thin curves show the paths corresponding to a disordered lattice for 16 realizations of disorder. 
V n is assumed to have a uniform distribution in the range (— k, k). The lower panels in (a) and (b) depict the distribution of 
the occupation probabilities \c n {t)\ 2 (in arbitrary units) at time t = 40/k for the ordered lattice and for one realization of the 
disordered lattice. The lower panel in (c) show the evolution of the normalized occupation amplitude a n (t) in a pseudo color 
map (modulus square of a n (t)). 

Fig-3 Wave packet evolution in a tight-binding lattice with two potential defects Vo at sites no and n\ (upper inset) for 
(a) the Hermitian lattice (h = 0), and (b) in the presence of an imaginary gauge field (h = 0.25). Parameter values are 
Vo/k = 1 and m — no = 20. Initial condition corresponds to the Gaussian wave packet c„(0) oc exp[—(n — 100) 2 /64 + inn/2]. 
The figures show the evolution of the normalized amplitude probabilities a n (t) (modulus of a n (t)) is a pseudo color map; the 
two vertical dashed lines in the maps show the position of the potential defects. While in the Hermitian case the two potential 
defects cause multiple wave packet reflection back and forth and echoes of the transmitted wave packet, multiple reflections 
and echoes are suppressed when the imaginary magnetic field is applied. 

Fig.4 (a) Numerically-computed quasi-energy spectrum (real and imaginary parts of the quasi energies E±) of the non- 
Hermitian driven Hamiltonian (7) in the absence of disorder for parameter values lo/k = 1, F/oj = 2, Gr/k = 4.7 and 
Gi/k = 4.26. The filled dotted curves correspond to forward propagating modes ( v g = R c(dE/dq) > 0) and positive imaginary 
quasi energy, whereas the open dotted curves correspond to the backward propagating waves. The dashed curve in the upper 
plot shows, for comparison, the quasi energy spectrum in the Hermitian case Gi = 0. (b) Numerically-computed evolution of 
the wave packet center of mass ( n(t )) in the semi-infinite Hermitian lattice corresponding to initial excitation of the edge site. 
The dashed curve shows the path in the absence of disorder, whereas the solid thin curves show the paths corresponding to a 
disordered lattice for a a few realizations of disorder. V„ is assumed to have a uniform distribution in the range (—k,k). (c) 
Same as (b), but for the non-Hermitian lattice. 

Fig.5 Optical realization of an imaginary gauge field in coupled microresonators, (a) Schematic of the CROW structure, 
comprising a sequence of main microrings indirectly coupled via auxiliary rings. The auxiliary rings provide light amplification 
(attenuation) in the upper (lower) half perimeter of the microring, (b) Schematic of indirect coupling of two main microrings 
via an auxiliary microring. 

Fig.6 (a) Schematic of a CROW laser made of N coupled microrings, (b) Distribution of resonance frequency detunings 
arising from disorder in a chain made of N = 60 microrings (behavior of U n /ft). (c),(d) Intensity distribution (in arbitrary 
units) of the lowest threshold lasing mode in the Hermitian (h = 0, panel c), and non-Hermitian (h = 0.05, panel d) cases. 

Fig.7 Numerically-computed wave packet evolution in (a) the Hatano-Nelson Hamiltonian [Eq.(4)], and (b) in the driven 
lattice Hamiltonian [Eq.(8)j in the absence (left panels) and in the presence (right panels) of disorder. The strength of disorder 
is as in Figs.2 and 4. The figures show the evolution of the normalized amplitude probabilities a n (t) (modulus of a n (t)) is a 
pseudo color map. Initial condition corresponds to a two-humped wave packet c n (0) oc n exp[— (n/10) 2 + iqon] with carrier 
Bloch wave number qo = 7r/2 (upper plots) and qo = 0 (lower plots). 
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Fig-3 
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